function [points, triangles] = triangulation(geometry, edge_length)

max_x = max(geometry(1,:));
max_y = max(geometry(2,:));
min_x = min(geometry(1,:));
min_y = min(geometry(2,:));
bound = [min_x max_x;min_y max_y];

[points,triangles]=distmesh2d(@dpoly, @huniform, edge_length, bound, geometry', geometry');
end